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ANALYTICAL STUDIES OF LUNI-SOLAR EFFECTS ON THE 
MOTION OF ARTIFICIAL SATELLITES 

By 

T. L. Felsentreger and J. P. Murphy 
Goddard Space Flight Center 


Abstract 

Theoretical and experimental studies indicate that long periodic 
lunar and solar forces produce perturbations in the orbital elements 
of some close earth satellites which are comparable to long period 
geopotential zonal harmonic effects. Analytical formulations for third 
body, solar radiation pressure, tidal, and precession and nutation effects 
have been successfully used to account for observable long period vari- 
ations in the orbital elements of such satellites as Relay 1, Relay 2, and 
Telstar 2. For example, in the case of Relay 1, a variation in the eccen- 
tricity with an amplitude of about .0007 and a period of about 400 days was 
successfully attributed to luni-solar forces using these formulations. 

Remaining variations in some orbital elements (after removal of the 
luni-solar, etc. effects) provide strong evidence for adjustment of geo- 
potential parameters (used in the orbit determinations) to more recent 
values. For example, use of the value -2.504 x 10 ~ 6 for J 3 would ac- 
count for the remaining variation in the eccentricity of Relay 2 (the value 
-2.285 x 10 ~ 6 was used in the orbit determination). 

Preliminary attempts to deduce seasonal variations in the geopotential 
have been promising. This study involves an analysis of the long period 
variations remaining in the orbital elements of satellites after removal of 
all zonal harmonic, luni-solar gravitation, solar radiation pressure, ti- 
dal, and precession and nutation effects. 

Discussion 


Introduction 


Most orbit theories for close earth satellites neglect lunar and 
solar forces with respect to the effects caused by the Earth's gravity 


field. However, an examination of the parameters in the long period 
parts of the disturbing functions indicates that some luni- solar forces 
give rise to perturbations which are comparable to long period zonal 
harmonic effects. Moreover, these perturbations are observable and 
predictable in the orbital element histories of some earth satellites over 
time spans covering hundreds of days. 

In order to demonstrate this, analytic expressions for lunar and 
solar gravitational (secular and long period), solar radiation pressure 
(long period), tidal (long period), and precession and nutation (long 
period) effects have been formulated (References 1, 2, and 3). These 
formulations were obtained from a straightforward integration of the 
variation equations for the Keplerian orbital elements (semi-major axis, 
eccentricity, inclination, mean anomaly, argument of perigee, and longi- 
tude of ascending node). The elements were referred to the equator of 
date and the equinox at an initial epoch; the integration was performed 
by assuming only the arguments of trigonometric terms to be functions 
of time. Resonances were treated as special cases. 

The satellite elements used in the study were obtained from a reduc- 
tion of tracking data using an orbit determination computer program 
based upon the 1959 artificial satellite orbit theory of Brouwer (Refer- 
ence 4). Hence, these elements are the so called ''mean/' or "double- 
primed," elements of Brouwer. The theory involves only the geopoten- 
tial zonal harmonics through J ; the following values were used: 
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J 2 = 1.08219 x 10“ 3 

J A - -2.123 

x 10“ 6 

J 3 = -2.285 x 10" 6 

J 5 = -2.32 

x 10- 7 . 


The study was carried out on a number of satellites. However, only 
some of the more interesting results for the satellites Relay 1, Relay 2, 
Telstar 2, and Nimbus 2 will be reported here. There was a particularly 
interesting resonance phenomenon in the case of Relay 2 for which a 
special analytic formulation was derived (Reference 5). Also, an ob- 
servable long period variation in the semi- major axis of Telstar 2 was 
found to be attributable to an indirect solar radiation pressure effect 
caused by passage of the satellite in and out of the Earth f s shadow 
(Reference 6). 
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Perturbations in Orbital Elements 


The secular and long period perturbations in the Keplerian orbital 
elements were obtained from integration of the following variation equa 
tions, which include coupling effects with zonal harmonics: 


de 

dt 


v 1 -e"^ 
na 2 e 



di 

dt 


cos i 

na 2 / 1 - e ^ sin i 



1 

na 2 /T - 


sin 


BR 

\ dd 




MD 




• • 

d4 _ _ _2_ BR _ 1 ~ e 2 2L + 4L 8e + f!£ Si , 
dt na Ba na 2 e 
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dt 

dfi 

dt 


COS 1 


BR /l - e * BR do> do> . . 

- + + — S e + — 8 x ,^(2) 


na 2 /l - e 2 sin i na 2 e '^ e 


1 


na 


BR d n e dft . 

— - — + - — 8e + — Si . 

2 /I - e 2 sin i < * 1 
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The integrations were performed by considering the ar guments of trigon- 
ometric terms to be linear functions of time - hence, resonances had to 
be treated as special cases. In. Equations (2), Se and Si are the perturba- 
tions obtained by integration of Equations (1), and t, i , and fi are given 
by (see References 1 and 4) 


4 = n < 1 - 


3J 2 a . 2 


4a 2 (1 - e 2 ) 3/2 


3 j 2 a 4 

(l-3cos 2 i) + - — e [10-25e 2 + 16v / l-e 2 
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n 


3J 2 « 2 , 
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[3(4 4 3e 2 ) - 18 (8 + 7e 2 ) cos 2 i + 7 (28 + 27e 2 )cos 4 i] , 


[4 - 9e 2 + 1 2/ 1 - e 2 - (40~5e 2 (3) 


15J..a 4 (2 + 3e 2 )cosi 

+ 36v / l~e 2 )cos 2 i] - / * 7 13 2 


32a 4 (1 - e 2 ) 4 


(3 - 7 cos J i ) 


Also, in Equations (1) and (2), the function "R" represents, in general, 
any of the secular and/or long period "disturbing" functions due to lunar 
and solar forces. An equation for the semi- major axis does not appear 
because only secular and long period perturbations are being considered 
here. Perturbations in e , i, i, «u, and fi arise from lunar and solar 
gravitation and solar radiation pressure; tidal forces and precession and 
nutation cause effects in i , ^ , and ^ . 

Comparison with Geopotential Zonal Harmonics Effects 


One can obtain an estimate of the magnitude of the secular and long 
period luni-solar perturbations relative to those due to the Earth's 
zonal harmonics by a study of the parameters appearing as factors in 
the disturbing functions. 

It can be seen from the 1959 paper of Brouwer (Reference 4) that 
the long- period part of the Hamiltonian F* p appears with the small fac- 
tors J 2 /a 5 , J 3 /a 4 , J 4 /a s , and J 5 /a 6 , and the secular parts F la and Fj a 
have J 2 /a 3 , J 2 /a 5 , and J 4 /a 5 as small factors (F* » although second 
order, gives rise to first-order long- period perturbations). The secular 
and long-period gravitational disturbing functions of the sun and moon 
have the factors m Q n 0 2 a 2 and m g n g 2 a 2 , respectively. The long-period 
solar radiation pressure disturbing function appears with -Fa as a mul- 
tiplier (for illustration, the value of A/m for Telstar 2, which is about 
0.074860427 cm 2 /gm, was chosen for use in this analysis). 
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Table l gives a listing of these parameters for increasing values of 
a, here the values of the zonal harmonic coefficients used w^rj those 

“ dopted by Goddard Computing Center (see Introduction). In addition 
the radius of the Earth has been chosen as the unit of distance, time is ’ 

radW VanBUard unit < 806 ' 832 sec.), and angular measure is in 

Assuming that divisors of approximately the same order arise in the 
tegration of both the zonal harmonic and luni-solar disturbing functions, 

Table 1 

Parameters in the Secular and Long- Period Disturbing Functions 

(Telstar 2) 


a 

(oartl 

radii 


J// a’ 

J,/a« 

-J./a 8 

J s /a 6 

1.00 

10.082x10 

11.70x10 

7 -22.85x10 

-21.23x10 7 

-232.0x10 

1 .05 

9.347 

9.173 

-18.81 

-16.64 

-173.1 

1.10 

8.129 

7.269 

-15.62 

-13.18 

-131.0 

1.15 

7.111 

5.8 20 

-13.06 

-10,56 

-100.3 

1.20 

6.262 

4.705 

-11.0? 

-8.535 

-77.72 

1.25 

J 5.540 

3.h:jo 

-9.360 

-6.957 

-60.82 

1 .30 

4,925 

3.153 

-8.000 

-5.718 

-48.06 

1.33 

4.398 

2.611 

-6.879 

-4.735 

-38.33 

1.40 

3.943 

2.176 

-5.948 

-3.947 

-30.81 

1.45 

3.549 

1.826 

-5.169 

-3.312 

-24.96 

1.5 

3.206 

1.542 

-4.513 

-2.796 

-20.37 

1.6 

2.642 

1.117 

-3.487 

-2.024 

-13.83 

1.7 

2.202 

.8244 

-2.730 

-1.495 

-9.611 

1.8 

1.855 

t 

.6196 

-2.177 

-1.123 

-6.822 

1.9 

1.577 

.4727 

-1.753 

-.8574 

-4.932 

2.0 

1.352 

.3657 

-1.428 

-.6635 

-3.625 

2.2 

1.016 

.2272 

-.9755 

-.4120 

-2.046 

2.4 | 

.7827 

.1470 

-.6888 

-.2666 

-1.214 

2.6 

.6156 

.09856 

-.5000 

-.1787 

-.7510 

2.8 

.4929 

.06801 

-.3718 

-.1231 

-.4814 

3.0 

.4007 

.04816 

-.2821 

-.08737 

-.3182 

3.5 

.2524 

.02229 

-.1523 

-.04042 

-.1262 

4.0 

.1691 

.01143 

-.08925 

-.02073 

-.056G5 

4.5 

.1187 

.006342 

-.05572 

-.00151 

-.02794 

5.0 

.08656 

.003746 

-.03656 

-.006794 

-.01485 


m t n< 2 a 2 


o.r>r»o7xio 

.6071 
.6663 
.72811 
.7930 
.8605 
.9307 
1.004 
1.079 
1.158 
1.239 
1.410 
1.592 
1.784 
1.988 
2.203 
2.605 
3.172 
3.723 
4.317 
4.956 
6.746 
8.811 
11.15 
13.77 


m , n 2 a 2 


-Fa 


2.579x10 
2.843 
3.121 
3.411 
3.714 
4.030 
4.359 
4.700 
5.055 
5.422 
5.803 
6.602 
7.453 
8.356 
9.310 
13.32 
12.48 
14.86 
17.43 
20.22 

23.21 
31.59 
41.26 

52.22 
64.48 


r " 3.537x10'* 
3,714 
3.891 
4.068 
4.244 
4.421 
4.598 
4.775 
4.952 
5.129 
5.306 
5.659 
6.013 
6.367 
6.720 
7.074 
7.781 
8.489 
9.196 
9.904 
10.61 
12.38 
14.15 
15.92 
17.69 
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we can see that at about a a 1.2 the long- period effects due to luni- 
solar gravitation are about the same as those arising from J 5 . Further, 
at a ~ 1.9 the luni-solar long-period effects are just as important as 
those due to any one of J 2 2 , J 3 , J 4 , or J 5 . Also, at a ~ 3, solar radiation 
pressure is as significant as the zonal harmonics of the Earth. For in- 
creasing values of a, of course, the luni-solar forces became progressively 
larger while the effects due to the Earth's zonal harmonics diminish. 

As far as the secular terms are concerned, the luni-solar effects 
became comparable to the zonal harmonic effects at about a % 5. There- 
fore, the integration becomes invalid at this point, and the formulas 
should not be used for satellites whose orbits have semimajor axes ex- 
ceeding this value. 

Analysis of Relay 1 and Telstar 2 


An examination of the mean orbital elements of Relay 1 and Telstar 
2 reveals quite substantial long-period variations remaining in the ec- 
centricity, inclination, argument of perigee, and longitude of ascending 
node. These elements are the "double-primed" variables of Brouwer 
(see Reference 4), so the effects due to the earth's zonal harmonics 
though J 5 have been accounted for. Since the two satellites are high 
enough so that air drag may be neglected, it was felt that these pertur- 
bations might be due to luni-solar forces. For these two satellites tidal 
and precession and nutation effects were found to be negligible compared 
to lunar and solar gravitational effects and solar radiation pressure. 

Data for Relay 1 and Telstar 2 spanning 630 and 525 days, respec- 
tively, were analyzed (Reference 1). Figures 1 to 4 give plots of the 
uncorrected and corrected values of eccentricity and inclination verjus 
time. It is clear that the major portions of the perturbations were indeed 
due to luni-solar forces. 

The argument of perigee and the longitude of ascending node were 
handled in a somewhat different fashion because of the presence of secu- 
lar effects. First, the existence of long-period variations in the double- 
primed elements was detected by examining the residuals from least- 
squares fits to a constant plus a secular term. Then the long-period 
luni-solar perturbations were subtracted from the double-primed ele- 
ments, additional least squares analyses of the same type were made, 
and the residuals were again examined. The least squares results are 
as follows: 


6 


INCLINATION (dag) 


« nvi< ■ 

0.28500 

0.28480 

£ 

5 0.28460 
z 

UJ 

U 0.28440 

UJ 

0.28420 
0 28400 • 
0 2B380 ( - 


• MEAN VALUE' * i-i 

* MEAN /AIUE5, LESS LONG -PERIOD LUNAR AND V'UR GRAVITATION/.! 
t f FECT 'j AND SOUR RADIATION rRfSSUSf it, I n 


• • 


* •• 


•.I 

..... 








40 80 


_t i__x_j.__.i_ 



t * >0 I doyi) 


NASA OSFC t AOS 

MISSION A TAAJCCIOAV ANAI.VSIS DIVISION 

MNANCH JU . OAU ..ItJxTl 

■ V /.^l.^., nOtNO J1»4 


Figure 1. Eccentricity of Relay 1. 


47.540 

47,530 

47.520 

47.510 

47.500 

47,490* 


• MEAN VALUES (!) 

« MEAN VALUES, LESS LONG - PERIOD LUNAR AND SOLAR GRAVITATIONAL EFFECTS 
AND SOLAR RADIATION PRESSURE (I, ) •, 


t- » 


... «»■’ 
* . I. : .» 


**»*■«>• «■ 


«>■ « «« * 




• ••• 


*** *e*i 


I I 1 — J i I I L— L_i 1 -j I 


. I — .1 i 1 — L—J L 


40 80 120 160 200 240 280 320 360 400 440 480 520 560 600 640 

... I Xu. \ 

' *0 V ®7* > MISSION A TAAMCtOflV ANALVtlt DIVISION 

SAANCM M2 OATS HUH 

•» i ma a m. **•«>? no vm — 


Figure 2. Inclination of Relay 1. 


7 



0 . 401 *) | — 

- • MIAN VALUES ( 

0.401,-3 - " MEAN VALIA:,, L ,$ LONG * PERIOD LUNA* AND SOU* GRAVITATIONAL 
IFFICTS AND SOU* RADIATION PRESSURE (•«)•• 

• •• 

0.40150 - . 


« - • wt 

, 1( , j b “-"-"“"""""■i"' 

■ ■ 1 * •_ 


• • 


0 40 80 120 160 200 240 210 320 360 400 440 410 520 560 600 

L , l UUWC TAM 

T 1 0 ( <nyi / MIMIOM A TAAMCTOAV ANM.VWA OIVIWON 

MAWCH _.m 0«r«... .Hl/M , 


Figur# 3. Eccentricity of Tolftor 2. 


• MEAN VALUES <l) 

■ MEAN VALUES, LESS LONG -PERIOD LUNAR AND GRAVITATIONAL EFFECTS 
AND SOLAR RADIATION PRESSURE ( l« ) , • 


, , «« >«> „i «** -- 


• • • » • 
• • 



120 160 200 240 280 320 360 400 440 400 520 560 600 

I - tn ( dav» ) UUTVC T1H 

•0 V I MMKW A TBAACTOAV ANALYtM MVMtON 

IAAWCH W OAU JUUn . 

■y »» »— w » n or wo. mi 


Figure 4. Inclination of Tel star 2. 


8 











RESIDUAL (dog) RESIDUAL (deg) 



l-l 0 (doyi) 


NASA-GSFC-TSDC 

MISSION * TRAJf CTOflV ANALYSIS OlVIStOM 

MIAIiTM OATi iLWI 

■ V ftl mn t ntw PL OT NO . » St 


Figure 8. Residuals from least squares fits to longitude of ascending node for Telstar 2. 





(1) Relay 1: 


= 3492704 + (122126 4/day) (t - t 0 ) 
w c = 8a; = 3492850 + (1221206/day) (t - t 0 ) 

R = 372035 - (1227953/day) (t - t 0 ) 

Q c = Q - SO = 362885 - (1227895/day) (t - t 0 ) 

(2) Telstar 2: 

" = 2972183 + (1221685/day) (t - t 0 ) 

— oj ~ 8ai = 2972242 + (1221716/day) (t - t Q ) 

n = 762503 - (1205 483/day) (t - t 0 ) 

n c = Q - 8fi =762446 - (1205 4 9 7/day) (t - t 0 ) 

Figures 5 to 8 indicate that the principal remaining variations in the argu 
ment of perigee and longitude of ascending node were caused by the sun 
and the moon. 

To illustrate the fact that long-period luni-solar effects can be as 
significant as first-order long- period zonal harmonic effects, let us ex- 
amine the perturbations in the eccentricity of Relay 1 and Testar 2. 

Tables 2 and 3 present the amplitudes of the principal long- period terms 
appearing in the formulas for the luni-solar and zonal harmonic effects, 
along with the corresponding arguments r nd periods (the formulas for 
the zonal harmonic terms were obtained from Reference 4). It is clear 
that the luni-solar effects are just as important as zonal harmonic 
perturbations. 

Analysis of Relay 2 


Similar results were obtained in the case of Relay 2. However, the 
situation was complicated by the existence of a near- resonant condition 
caused by a longitude of perigee (o> + R ) which was practically constant 
for the time interval studied (about 654 days). The period of the longitude 
of perigee was found to be on the order of 550 years. This condition 
caused an almost secular increase in the eccentricity and, to a much lesser 
extent, in the inclination. 
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Table 2 

Principal Long- Period Terms in Zonal Harmonic and Luni-Solar 
Perturbations in Eccentricity of Telstar 2. 


Source 

Amplitude 

Argument 

Period (days) 

(Jj/J j) 

373 x 10 « 

■* 

296 

(J4/J2) 

-60 

2-* 

148 

(1.) 

17 

2a, 

148 

(Jj/Jj) 

5 

U 

296 

(J5/J2) 

-1 

3u 

99 

Lunar gravity 

-220 

2{% - *-«) 

1121 

Lunar gravity 

192 

Q e - 2a - t) 

261 

Lunar gravity 

-95 

2a 

148 

Lunar gravity 

33 

2(A. fi “ a - >i) 

14 

Lunar gravity 

11 

2K C - !i c ' 2a - £1 

14 

Lunar gravity 

-11 

n c * 2a - n 

103 

Solar gravity 

245 

2(\, -a-n) 

219 

Solar gravity 

221 

2A id * " 2a - 

607 

Solar gravity 

-106 

• 2 (.\-a.-n) 

1111 

Solar gravity 

90 


261 

Solar gravity 

-42 

2a 

148 

Solar gravity 

-35 


781 

Solar gravity 

8 

' u, o ‘ 2u; " 0 

917 

Solar gravity 

-6 

♦ 2a - fl 

103 

Solar radiation pressure 

42 

OJ + Q ~ k @ 

437 

Solar radiation pressure 

-24 

w ■ ^0 

1562 
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Table 3 

Principal Long- Period Terms in Zonal Harmonic and Luni- Solar 
Perturbations in Eccentricity of Relay 1. 


Source 

Amplitude 

Argument 

Period (days) 


520 x 10 " 6 

CO 

298 

{]>/!,) 

-71 

2io 

149 

(h) 

18 

2oj 

149 

(h/h) 

7 

O) 

298 

(J./J.) 

7 

3u> 

99 

Lunar gravity 

304 

2(n c -a-fl) 

2584 

Lunar gravity 

147 

n c - 2 oj - n 

315 

Lunar gravity 

-68 

2a, 

149. 

Lunar gravity 

19 

2{K c -*-n) 

14 

Lunar gravity 

9 

2 \ c - n c - 2 ^ - n 

14 

Lunar gravity 

-7 

n c + 2 a. - n 

97 

Solar gravity 

149 

2 (n Q -c-n) 

2691 

Solar gravity 

109 

2(\ s - a.' - n) 

171 

Solar gravity 

101 

2 k„ - n e - 2 a, - n 

436 

Solar gravity 

70 

- 2w - n ' 

314 

Solar gravity 

-30 

2a 

149 

Solar gravity 

-30 


795 

Solar gravity 

12 

~ n 

2256 

Solar radiation pressure 

31 

a; + fi - 

342 

Solar radiation pressure 

-26 


1590 











The luni-solar gravitational and solar radiation pressure effects 
were computed using the formulas in Reference 1. However, the pre- 
caution was taken to exclude the near- resonant terms with arguments 
2 (U - - 12) and 2(i '1% - - 0) from the initial computations. The por- 

tion of the luni-solar gravitational disturbing function which includes the 
two near- resonant terms is 


Fr ~ 6 ? * 2q2 ^ + C ° S i ' )2 ^ n ® m& Sin2 COS ” a ' “ n > 


+ n 2 sin 2 cos 2(0. - cu - f2)] . 


The method of integration employed in Reference 1 cannot be directly 
applied to the "moon" term in its present form inasmuch as the motions 
of i,£ , f) ?f and oo + 0 are all commensurate; in addition, since i, e and fig 
have been referred to the Earth's equatorial plane, their motions are not 
sufficiently constant over the time interval under study to allow accurate 
computation of the perturbations. 

However, using the ecliptic as a reference plane the "moon” term 
can be rewritten in terms of angular variables having motions which are 
essentially constant, and the integration can be carried out (details may 
be found in Reference 5). Figure 9 illustrates the effect of removing the 
luni-solar perturbations, including the "near- resonant” perturbation (^e) R , 
from the "mean” values of eccentricity. 

There still remained a variation in the eccentricity having a period 
approximately that of the argument of perigee. A harmonic analysis us- 
ing the corrected values of the argument of perigee (a r5 ) was performed, 
the results of which are shown in Figure 10. A further analysis indicated 
that this variation could be accounted for by use of the value -2.504 x 10~ 6 
for the J 3 zonal harmonic, instead of the value used in the orbit deter- 
mination (see Introduction). This "adjusted” value for j 3 is much more 
in line with recent determinations. 
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Figure 9. Eccentricity of Relay 2. 






Shadow Effects on Telstar 2 


An indirect long period solar radiation pressure effect on the orbit 
of a satellite is caused by the passage of the satellite in and out of the 
earth's shadow. This effect is best observed in the semi-major axis, 
which undergoes no long period motion from any other source. 

The Telstar 2 satellite provides a good subject for a study of this 
effect. A long- period variation with approximate amplitude of 0.0000140 
earth ra<ti and a period of about a year was observed in the semi-major 
axis. Analytic formulations derived using a procedure described in 
Reference 6 successfully accounted for this variation, illustrated in 
Figure 11. 

Seasonal Variations in the Geopotential 


A seasonal variation in the second zonal harmonic in the geopoten- 
tial has been observed (Reference 8). In his paper, Dr. Kozai sought but 
did not find a variation in the third zonal harmonic. An analysis is under- 
way to use the moderate to highly inclined orbits of the close nearly 
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Figure 11. Semi-Major Axis of Telstar 2 vs. Time 
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circular Tiros 8, Alouette 1, and Nimbus 2 satellites to study these ef- 
fects. Orbits such as these undergo very substantial periodic perturba- 
tions due to the third harmonic in eccentricity and argument of perigee 
(see Reference 7), Any residual periodicities after all the luni-soJar 
gravitation, precession and nutation, radiation pressure, high degree 
zonals, and solar and lunar tidal effects are removed may reveal the 
effects of the time dependent geopotential. 


Symbols 

A 

a 


e 



4 

l 

m 

n 


Effective presentation area of satellite 
Semi-major axis of satellited orbit 
Mean equatorial radius of the Earth 
Eccentricity of satellite's orbit 

Eccentricity corrected for long-period lunar and solar 
effects 

Solar radiation pressure force constant 
Resonant disturbing function 

Inclination of satellite's orbital plane to Earth's equa- 
torial plane 

Inclination corrected for long-period lunar and solar 
effects 

Zonal harmonic coefficients in Earth's gravitational 
potential 

Mean anomaly of satellite 
Mean motion of mean anomaly 
Mass of satellite 
Mean motion of satellite 
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R Disturbing function 

t Time 

a - Argument of perigee of satellited orbit 

• 

^ Mean motion of argument of perigee 

n Longitude of ascending node of satellited orbit 

H Mean motion of longitude of ascending node 

Mean longitude of sun or moon 

The symbols e 0 , i g , etc. refer to the corresponding quantities for the 
sun (©) or the moon («). 
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